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Passive turbulent fiamelet propagation 

By Wm. T. Ashurst, 1 G. R. Ruetsch 2 AND T. S. Lund 2 

We analyze results of a premixed constant density flame propagating in three- 
dimensional turbulence, where a flame model developed by Kerstein et al (1988) 
has been used. Simulations with constant and evolving velocity fields are used, 
where peculiar results were obtain from the constant velocity field runs. Data from 
the evolving flow runs with various flame speeds are used to determine two-point 
correlations of the fluctuating scalar field and implications for fiamelet modeling axe 
discussed. 


1. Introduction 

For some applications, it is useful to consider the premixed flame structure as 
unchanged by turbulence and to cast the problem as a fiamelet moving through 
the flow with a known density jump (Linan & Williams, 1993). A further assump- 
tion is to ignore the density change and consider passive fiamelet propagation. We 
find that these constant-density flames do reveal possible mechanisms of turbulent 
flame propagation (Ashurst, 1994). Therefore, the fiamelet assumption has been 
exploited in order to decouple the complexity of chemistry from that of turbulence. 
To investigate flame stability, Markstein wrote a flame evolution equation in 1964, 
but lacking a computer, he was restricted to a stability analysis. Markstein’s no- 
tation was / for a flame with volume expansion; while a passive formulation has 
acquired the letter G\ we distinguish between these fiamelet models by using / when 
volume expansion is considered and G when it is not. The special initial condition 
formulated by Kerstein et al (1988) is G = x for the evolution equation 

^ + u-VG = 5 l |VG| (1) 

where the right side describes Huygens’ propagation with Sl as the burning velocity. 
This nonlinear term makes the initial condition very special: it provides a connection 
between flame area and the gradient of G. The scalar gradient corresponds to the 
flame surface density, that is, flame area per unit volume. In this passive formulation 
any level surface represents a flame and spatial derivatives of G define the geometry 
of that particular surface, thus the actual magnitude of G is not important. This 
formulation allows a numerical simulation which follows Damkohler’s idea: creation 
of flame area by turbulent motion causes an increased consumption rate. Now, 
with the passage of fifty years, the computer can easily determine passive flame 
area within a turbulence simulation and, when the motion is forced, a statistically 
steady propagation is achieved. 
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FIGURE 1. Autocorrelations ;jf?n (top) and 1 R 22 (bottom) for u f = Sz,/2. 
:i = 1, :i — 2, and :i = 3. 

2. Summer objective 

The initial suggestion for the Summer Program was motivated by the previous 
work of Lund & Rogers (1994) which examined the eigenvalues of the strain-rate 
tensor in forced, homogeneous, isotropic turbulence. Their velocity fields, on meshes 
up to 256 2 3 , could be used as a frozen flow in which the propagation of the scalar 
G would be done, and thereby, obtain the G statistics at a smaller cost. This 
suggestion did not work. The frozen flow created very large distortions in the G 
level surfaces, so large that a single surface was connected across the periodic system 
length. This occurred because the amplitude of flame distortion in a frozen flow 
is proportional to the rms velocity of the velocity mode (Ashurst ei a/., 1988) and 
the lowest modes have the largest energy content in these forced simulations. Since 
any level surface represents a flame in this passive model, the large distortions 
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FIGURE 2. Autocorrelations iR\\ (top) and 1 R 22 (bottom) for u 1 = 2 Si. The line 
types follow the convention in Fig. 1. 


correspond to flames interacting with their own periodic image and this is not 
desirable. The distortion amplitude is reduced when the flow has a time variation, 
and so evolution of both the velocity field and the scalar G field has been done with 
32 3 systems. 

Previous G field simulations have been done with finite-difference techniques 
where the truncation error combined with grid-resolution suppresses large gradi- 
ents in the G field. These large gradients correspond to a cusp in a G surface, and 
the cusps are formed by the Huygens’ propagation mechanism. Cusp formation and 
the pseudo-spectral solution technique are not compatible, and so a hyperdiffusivity 
term was added to the dynamical equation of the form £)V 4 G with D = \v. 





154 


Wm. T. Ashurst , G. R. Ruetsch & T. S. Lund 



r 



r 


FIGURE 3. Cross-correlations iR \2 for u ' = Sl/ 2 (top) and u' = 2 Si (bottom). 
The line types follow the convention in Fig. 1. 

3. Isotropic propagation 

Are the scalar fluctuations of G isotropic? Peters (1992) considers an ensemble 
of flamelets travelling from all directions into a spatial domain, and so within this 
domain the fluctuations are assumed to be isotropic. Thus, the well- developed 
analysis of isotropic, homogeneous turbulence may be applied to the geometry of 
these flamelets. The geometry is related to the scalar fluctuations as a consequence 
of the initial condition that G = a\ The fluctuations correspond to the negative of 
the flame displacement from the transverse plane associated with that particular 
flame: use g f = G — x, find by adding a constant so that G = 0, we obtain g* as the 
negative displacement from the mean flame location where G = 0. The additive 
constant can be done for any plane normal to the mean propagation direction, and 
thereby we relate g ' values to flame displacement from a transverse plane. 
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FIGURE 4. Passive front propagation through the swirling flow created by a 
Gaussian distribution of vorticity, the maximum swirl velocity is 0.8 Si and the 
spatial unit is 4 r) with a time interval between flame images of The filled 

circles are points which are moved in the local flame normal direction, and the 
change in their spacing with time indicates the stretching of the flame surface. 

To examine this issue with the numerical simulations, two scalar fields were com- 
puted within the domain: G\ with mean gradient in the x direction and G 2 with a 
gradient in the y direction. For analysis of the scalar fluctuations, the mean gradi- 
ent is removed (<?i = G\ - x + x Q and p 2 = G 2 — y + Vo) and the values of x Q and 
y 0 are determined so that g\ and <72 have a zero mean volume average. Division by 
their rms values provides autocorrelations with value of unity for zero separation. 
We can define a general two-point correlation as: 

iRjk{r) = gj(x)gk(x + re<) 

where autocorrelations refer to instances where j = k and cross-correlations other- 
wise. Autocorrelations for simulations with u ' = Sl/ 2 and u' = 2 Sl are shown in 
Figs. 1 and 2. This data in these figures are taken from a single realization of the 
scalar field. Due to the limited sampling the correlations are not fully converged, 
but do exhibit dominant trends which we now discuss. 

When u 9 — Sl/ 2, the autocorrelation is larger in the direction of the imposed 
mean gradient (for either (/i or ^ 2 )- Therefore, the Huygens’ propagation domi- 
nates over the distortions created by the velocity fluctuations and the passive scalar 
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fluctuations do have a directional bias. However, when u* = 2Sl, then the autocor- 
relations no longer show as strong a directional dependence. 

The cross correlations indicate a negative value for zero separation and an asym- 
metric nature with respect to direction of the separation, see Fig. 3. We now 
speculate on why this is so. In Fig. 4, we present flame shapes created by passive 
propagation through the swirling flow created by a Gaussian distribution of vortic- 
ity. The propagation is from right to left, corresponding to a scalar G field with 
mean gradient increasing towards positive x. To the right, before the flame interacts 
with the vortex, the straight flame implies that g* ~ 0 along that flame surface. On 
the left side, after the interaction, the flame is distorted so that above the vortex 
center g f > 0 and below the center g 9 < 0. Hence, the swirling flow combined with 
Huygens’ propagation creates asymmetric flame shapes. Taking Fig. 4 and rotat- 
ing it a quarter of a turn counter-clockwise gives the solution for the scalar field 
G 2 which has its mean gradient towards +y. Now the cross- correlation between 
Gi and G 2 > using g[ and g* 2 , has a significant contribution in only one quadrant, 
the lower-left quadrant in Fig. 4, where g[ < 0 and g' 2 > 0 producing a negative 
cross- correlation for zero separation. 

4. Future work 

Visualization of the intense vorticity and its relationship to flame area will be one 
aspect of future work. In previous work, Ruetsch & Maxey (1992) found that intense 
passive scalar gradients occur between vortical regions; will the same structure occur 
in a propagating, passive scalar field? 
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